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Identity for the second functional derivative of the density functional Hartree plus 

exchange-correlation functional. 
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It is shown that the second functional derivative of the density functional Hartree plus exchange- 
correlation functional satisfies 



= constant. 



Sp N (r') 8p N (r) 



Pn (r) and p 1 N _ 1 (r) are TV-electron and (N — l)-electron densities determined from the same adia- 
batically scaled Hamiltonian of the interacting electron system with 7 the scaling parameter of the 
^ . electron-electron interaction strength. 



INTRODUCTION 



The Kohn-Sham (KS) formulation of Density Functional Theory (DFT) has become the de facto tool for the 
calculation of electronic structure of molecules and solids. In all practical applications of DFT, however, approxima- 
tions to the exact functionals have to be made Properties of density functionals that give an indication of the 
internal structure of the functionals can give insight and help with thee development of accurate approximations to 
the exact functionals 

There is little known about the properties of the second functional derivatives of important functionals. From 
structural stability arguments jbM12l ]. for example, the second derivative of the density functional F 1 [p] , the sum 
of the kinetic and interaction energy functionals, is positive definite. Similarly for the non-interacting Kohn-Sham 
, system, the second derivative of the non-interacting kinetic energy functional is positive definite. In this paper it is 
shown that the integral of the product of the charge density difference of the N and (N — 1) electron densities of 
the same Hamiltonian and the second functional derivative of the Hartree plus exchange-correlation functional of the 
1^1 ' iV-electron density is equal to a constant: 

>: /^(^(O-A-W) jJZUgL -n. tan,. (1) 

Here pn and p' t N _ 1 are the ground state charge densities of an interacting TV and (N - 1) electron system of the 
same Hamiltonian with multiplicative external potential Ug Xt ([/O/v])- The potential Ug Xt ([/Cat]) is constructed to keep 
the charge density of the N electron system independent of the coupling strength parameter 7 [IB - [l6| that scales the 
electron-electron interaction strength. At 7 = 1 full strength Coulomb interaction between electrons is included and 
the external potential v^ xt ([pn]) is the external potential of the fully interacting system, while 7 = corresponds to 
the non-interacting Kohn-Sham potential. 

Equation (JTJ) gives new insight into the internal structure of density functionals and can also be used to test 
approximate density functionals. The simplest test can be performed at full coupling strength, 7=1. Two self- 
consistent calculations are required to determine pn and p]v-i an d then |T]) can be evaluated. Ideally an approximate 
$J] ■ exchange-correlation functional will satisfy this expression. 



II. PROOF 



In the adiabatic connection approach [13Hla | of the constrained minimization formulation of density functional 
theory @, 0, [l(| E3 the Hamiltonian H 1 for a system of N electrons is given by 

IT'=f + <yV ee +i)l tmt \p N ]. (2) 
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Atomic units, h = e 



1 are used throughout. T is the kinetic energy operator, 



T 



1 N 



(3) 



and ~/V ee is a scaled electron-electron interaction, 
The the external potential 



y cxt \PN 



N 



2<j 



(4) 



JV 



[pw] =X) u e X t([Piv];ri). 



(•5) 



is constructed to keep the charge density fixed at pjv (r) , the ground state charge density of the fully interacting 
system (7 = 1), for all values of the coupling constant 7. The external potential has the form [l7l [l8j 



«ext(biv];r) = (l-j)v hx ([p N ];r) 



«c([pjv];r) -^([pJvls^+^xtCbiv];^, 



(6) 



where Wg Xt ([piv] i r ) = w cxt (r) is the external potential at full coupling strength, 7 = 1, and Wo Xt ([p7v] ;r) is the non- 
interacting Kohn-Sham potential. The exchange plus Hartree potential [TTI. fl2j| vhx{[pn]] r) is independent of 7, while 
the correlation potential vJ{[pn\] r) depends in the scaling parameter 7. The chemical potential 



M = ^Wxt[p])-^-iWxtW) 



(7) 



depends on the asymptotic decay of the charge density [ll|, EH, Ell [20] and hence is independent of the coupling 
constant 7 [U [23]. In Eq. Q is the groundstate energy of the (JV — l)-clectron system with the same 

single-particle external potential Ug Xt ([pat] ; r) as the JV-electron system: 



H 1 

M 



*\ ) = E 1 



M 



Pi 



M 



^X/.oxt [pjv] = ^xt ([pjv] ; r 

The energy functional F 1 [p] is defined as [Tfj| - [l2l | 



(8) 



i=i 



= min ( \I/ 



T + -yV e 
f + 1 V c 



p 



(9) 



where according to the Levy constrained minimization formulation [Toj . the wavefunction 1^2) yields the density 



and minimizes ( W 



T + jVee ^/j . For w-representable densities [1CH12| |. the densities are derived from the groundstate 

eigenfunctions of the Hamiltonian in ([5]). Note that by construction of Ug Xt ([pn] ', r) , Eq. ^ pn — p\q is independent 
of 7, but the groundstate density of the (JV — l)-electron system p 1 N _ 1 is expected to be a function of 7. F^lp] is 
usually decomposed as [HI E2 

(10) 



F*\p]=T° [p]+ 7 S te W+^[p] 
The correlation energy EJ [p] is defined as [13] 



T + jV e( 



* 7 

(0 



If 1 



T + 7V* 



p 



(11) 



where |^°) is the Kohn-Sham independent particle groundstate wavefunction that yields the same 
interacting system at coupling strength 7. [p] is the sum of the Hartree and exchange energy 



and the kinetic energy functional T° [p] is given by 
The full kinetic energy 

T>\p] -- 



T 



T° [p] + T2 \p\ , 



with the correlation part of the kinetic energy defined as 

t 



T 



Assuming that f 17 [p] is defined for non-integer electrons [III 13 Hil , at the solution point 

+ «oxt(W; r ) =v 



5F^[p] , 



*p(r) 

where /i is the chemical potential. Note that by definition of F 1 [p] 

K {v2 xt [Pn]) = F~<[p N } + [ d 3 rp N (r) v^ t Qp N ] ; r) 



From and (fT7)) . 



Since 



^-1 ("«t = F^IPn-i] + / d ( r ) «cxt (bjv] ; r) • 



F"<[p N ] - F^lpjf^] =fi- d 3 r (p N (r) - p 7 r _ 1 (r)) w 7 xt ([pjv] ;r) 



tfV (r) - p 7 r _ 1 (r)) = 1, 



it follows from JIB]), (IS) and O that 



F 7 [pjv] _ Fi\pl_ x ] = / d 3 r (p N (r) - p 7 ^ (r)) 



6F~< [ PN ] 
5pn (r) 



In a recent paper [24{ the author showed that 

The proof of this equation is based on the virial theorem [llj of the interacting system. 
The charge density pn is independent of 7 by construction, therefore from (|10|) and ([12 

^-F^\p N } = E hx [p N } + ^-EZ{p N ]. 
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The (N — l)-electron density pj ! _ 1 is a function of 7, hence it follows from (fT0|) and (fT2|) that 

dpl_ t (r) <5r° [^_J 



+ - / d\ d -^^v hx ([pl^] ;r) + ^EJ[pl_ x ], (23) 



where u te ([p] ; r) = S ^ r \ p] ■ It follows from the definition of E2[p N ], Eq. CTTJ that [HI 

|^c 7 M - ^ - T2\p N ]) ■ (24) 

The derivative -g^EJ [p]q-i\ is l ess simple since pjj_ 1 ^ s a function of 7. From (fTT|) it is clear that -BJ [Pjv-i] depends on 
7 as the scaling parameter of the Coulomb interaction strength and via the dependence of the wavefunction 
on 7. It can be shown that[25[ 



Pn-i 



E2[pI_ x ]-t?IpI_ x ] 



7 



+ / dV ^^ 1 (r) (v2 ([pn] ; r) + jv hx ([p N ] ; r) - 7 v hx ([p^_J ; r)) (25) 
From (EH), (EH), (HHI), and 422} it follows that 

eLApn] - T2\p N ] - [eim^] - T2 [pl_,]) 

d\ ( PN (r) - pl_ x (r)) (El xc [p N ] - T2[p N }) 

^%lMff^l-f^V (26) 
9 7 \ fyjv-i W fyjv (r) J 

Here ^LcM = lE hx [p] + EJ[p}. From © and (TO]) 

7K 7 e [p] = sL c [piv]-r c 7 [p], (27) 

and therefore from (f2"6"|) and (j2~lj) 

A %w(gM.y-j). t (28) 

9 7 ^ (5pjv (r) Spjf^ (r) y 
By construction, from (|16p and (jTUJ) 

ST° [ PN ] 
Spw (r) 



+ u Lc([/'A f ]; r ) + ^xt(M;i') = M;v 



^p1_i 0) 



+ «L c (K-i]; r )+«c X t(M;r) = ^-i ( 29 ) 



where v^ xc ([p] ; r) = -npFp and hence 



(5^ (r) 5p^_ 2 (r) 

w L c ([pIt-i] ; r ) - u Lc(M; r ) + ^-p?v_ 1 . (30) 
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The dependence of p N _ x (r) on 7 is through the 7 dependence of the potential v2 xt ([pjf] ; r) = u 7 (r) . Using (|30| in 
k71> and since 



#7 



(iV - 1) = J £ 



dPN-i ( r ) 
c?7 



= 0. 



it follows 



<*v/ rf 3 r(,L c (W;r)-.Lc([^- 1 ] ;')) 



N 



<9u 7 (r') 
_i ^7 



(31) 



(32) 



The functional derivative with respect to v 1 is taken at constant particle number since pTv-i ano - Pn are derived from 
the groundstates of the Hamiltonian in Eq. © . 
The author recently showed that (2(| 



= J d\ ( PN (r) - pl_ x (r)) ^A-^ c[pjv ]. 
Take the functional derivative of (|33|) with respect to (r) and use he functional chain rule: 

/ ^r'(^c(M;r)-^c([P 7 W -i] !')) %^ 



(33) 



N-l 



= d 



'3 Jl I jZj 



[p N (r ) - p N _, (r )) — 



5p N (r") 5vt (r) 

^ v 7/ r "\ _____ 

Multiply both sides by — £ } and integrate over r . It follows from ([Ml) and (1571) that 



(34) 



where 



'V Ur"{ PN (r>)-pl_ 1 (T')) 



(JpiV (r) 
SpN (r) 



^ 7 (r") 
^7 



£«ext([PJv];r') 



(35) 



(36) 



JV 



is the density response function of the TV-electron system. From stability considerations Xn ( r > r ') i s negative semi- 
definite and has one zero eigenvalue which corresponds to the invariance of the density when the potential is changed 
by a constant [HI, HH. However, from © 



d 



— ^ — = -v hx ([p N ];r) - — v2([p N ];r) ^0, 
Hence, from the properties of Xn ( r > r ') Eq. (ESI) implies the main result of this paper, 

dV ( PN (rO-p^Jr')) 



5p N (r) 



= constant 



(37) 



(38) 



since 



Sp N (r) Sp N (r')Sp N (r) 



III. DISCUSSION AND SUMMARY 



In the derivation of Eq. (|3"B"|) use was made of the groundstate wavefunctions of the same Hamiltonian, for example 
in the derivation of (|25|) . Accordingly the relation (|38t has been proven for w-representable pure state densities [TTI.[l2l| 
only. Whether it is valid for non-u-representable densities and non-integer densities remains an open question which 
is under investigation. Throughout the assumption was made that all functional derivatives are well behaved and this 
implies that the functionals are defined for non-integer particle numbers as discussed in [TlL [l2l [23j . 
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In summary, an integral expression was derived for the second derivative of the pure state density functional Hartree 
plus exchange-correlation energy which provides new insight into the structure of density functionals. It also places 
constraints on potential approximate exchange-correlation functionals and can be used to test existing and potential 
approximate functionals. 
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